Single hole spectral function and spin-charge separation in the t — J model 
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Worm algorithm Monte Carlo simulations of the hole Green function with subsequent spectral 
analysis were performed for 0.1 < J/t < 0.4 on lattices with up to LxL = 32x32 sites at temperature 
as low as T = J/40, and present, apparently, the hole spectral function in the thermodynamic limit. 
Spectral analysis reveals a 5-function-sharp quasiparticle peak at the lower edge of the spectrum 
which is incompatible with the power-law singularity and thus rules out the possibility of spin- 
charge separation in this parameter range. Spectral continuum features two peaks separated by a 
gap ~ 4 -f 5 i. 
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For almost four decades the problem of hole dynamics 
in magnetic systems has attracted constant interest with 
applications ranging from properties of chaxge carriers 
in magnetic semiconductors and insulatoraUD to vacan- 
cies in solid ^HeH. The research in this area exploded 
with the discovery of high temperature superconductors 
in cuprates, where superconductivity appears upon light 
doping of AFM insulators. Despite an enormous theoret- 
ical effort over the years and quite a variety of treatments 
(for reviews, see, e.g., Ref. ^ a complete solution of this 
inherently strong-coupling problem still does not exist, 
especially in the most interesting region of i > J, where 
J is the exchange coupling constant and t is the hopping 
matrix element in the t — J Hamiltonian 

H^-tY^ ctcjs + jYl (s» - ■ (1) 

<ij>,s <ij> ^ ^ 

Here Cja is projected (to avoid double occupancy) 
fermion annihilation operator, = cj^ Cis 7^ 2 is 

the occupation number, S; = J2ss' "^L Css'Cis' is spin-1/2 
operator, and < ij > denote nearest neighbor sites of the 
2D square lattice. 

The central problem in the hole dynamics is whether 
or not its spin and charge degrees of freedom separate. 
The standard way to answer this question is to study 
the spectral function Ap{u!) — —TT~^Im Gp{uj), where 
G'p(ijj) is the hole Green function. If there is an ele- 
mentary excitation associated with the hole, the spectral 
function is supposed to feature a peak at the lower edge of 
the spectrum. What is crucial, however, is notithe pres- 
ence of the peak itself, but its functional formB. Within, 
the self-consistent Born approximation scheme (SCBA)III 
one finds finite overlap between the bare hole and low- 
energy quasiparticle states, which means that the peak 
is ^-functional and the hole is described as coherently 
propagating spin-polaron in the nearly ordered antiferro- 
magnetic (AFM) background (with vanishing scattering 
at low temperature due to small density of spin waves)Ll. 

In contrast to that, various resonating- valence-bond 
(RVB) descriptions and Anderson's general arguments 



about breakdown of the Fermi-liquid picture in the 
system with no-double-occupancy constraint (see, e.g., 
Ref. ^) strongly suggest that power-law singularity, which 
is indicative of spin-change separation scenario, might 
be the case (there is even a claim that the quasiparti- 
cle weight Z should be rigorously zeroa). To make the 
issue more confusing,. RjjgJe resolved photoemission spec- 
troscopy experimentstjO show very broad maximum in 
A{u}) which can be considered both as the quasiparticle 
peak with anomalously large broadening oii-as the evi- 
dence for composite nature of quasiparticlesEJ. 

The quasiparticle pictuE£_wp^ supported by exact cal- 
culations on small clustersQilj'Lj, but system sizes (up to 
32 sites) were too small to perform finite-size scaling. 
Variational calculations. Green function Monte Carlo 
and density matrix renormalization group studies were 
mostly concerned with the dispersion law Ck (lowest en- 
ergy in a given momentum sector). Large scale sim- 
ulations of the imaginary time Green function Gk{T) 
were performed recently using a combination of the loop- 
cluster Monte Carlo method for the AFM state aiid 
hole evolution in the fixed space-time spin backgroundEj. 
This method works only for relatively large exchange 
J > 0.6 t, since magnetic background is simulated with- 
out the hole and polaron-type distortions have to be ac- 
counted for as quantum fluctuations before the hole is 
introduced. For J/t < 0.6 the error bars in Gu{t) are 
too large for reliable spectral analysis (see below). 

To summarize, we still lack evidence that for small J 
the quasiparticle weight remains finite in the thermody- 
namic limit and the lowest peak has nothing to do with 
the power-law singularity. We thus find it important to 
rule out the possibility that t- J model may explain the 
data of Refs. |lO|,^ (as suggested by Ref. 0)- so that exten- 
sions of the model such as t' and t" termal3 or frustrating 
exchange couplings are proven necessary. 

Speaking classically, moving hole breaks AFM bonds 
and thua||-i±s energy increases linearly with the travel 
distanceBEZl (this consideration, or the string-potential 
picture, is most appropriate for the t—J^ model). It is be- 
lieved that the ground energy scaling £^ko J^^^ [where 
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ko = (77/2, ■7r/2)] and excitaticm-^pfictpum are described 
by the string-potential pictureQu'llSEjIlj, and transverse 
spin fluctuations do not "erase" strings completely. In 
the limit of small J the theory predicts that several res- 
onances in A{lu) have to be seen with the peak posi- 
tions being strictly related to the eigenvalue properties 
of Airy functions. Early exact diagonalization studies 
on clusters 4 x 4a attributed two peaks above the ground 
state tOjSlxijig resonances, however later studies on larger 
clustersoE3 were not able to detect the second resonance, 
and the spectral function was showing strong size depen- 
dence. What happens at small J in the thermodynamic 
limit remains an open question. 

In this letter we present results for G'k(T) and Ak(w) 
obtained from Monte Carlo simulations on systems with 
16 X 16 and 32 x 32 sites and at temperatures as low 
as T / J — 0.025 (for the largest system size) using 
continuous- time Worm algorithmcJ in combination with 
the recently developed spectral analysis which is. capa- 
ble of resolving infinitely sharp features in The 
method itself is free from any systematic errors, and 
we were unable to detect finite-size corrections in our 
data; thus, we believe, our results describe correctly the 
thermodynamic limit. In the parameter range studied 
0.1 < J/t < 0.4, the lowest peak in (w) is a 5-function 
within the resolution limit of order of 0.01-^0.03 t, which 
means that our quasiparticle is the spin-polaron. For the 
excitation spectrum we observe two well-separated peaks 
for all values of J. The ground state energy scaling does 
follow the J^/'^ law predicted by the string-potential the- 
ory. Although we were unable to resolve individual string 
resonances, we believe that their combined effect is seen 
as the first peak in the spectral continuum since its posi- 
tion also scales as J^/^. The high-energy peak is roughly 
at a constant distance ~ 5< from the ground state. 

Worm algorithm is based on the idea that world-line 
configurations of spins and the hole are updated through 
the space-time motion of the creation and annihilation 
operators. In Fig. 1 and Fig. 2 we show the typical config- 
uration of the Heisenberg AFM with the hole, and Monte 
Carlo updates which apply to it. Physical configurations 
contributing to the hole Green function are those which 
have no spin end points (denoted by filled circles). 

Since the world-line representation is based on the ex- 
pansion of the statistical evolution operator e~^^^ in 
powers of t and J it suffers from the sign problem which 
first appears in order t^J^ (see Fig. 3). It is worth noting 
that if not for the sign problem, spin-charge separation 
can be ruled out by the analysis of world-line configura- 
tions. Let the hole be created by Cii operator. If spin 
charge separation does take place, one should see, fol- 
lowing the evolution of the system configuration in imag- 
inary time, an extra spin density leaving the hole creation 
site and going into the bulk, i.e., the world- line density 
far from the hole should increase. At T = the AFM 
ground state is ordered (as opposite to the spin liquid 
state) and the minimal possible change in the world-line 
density is equivalent to having exactly one extra world 




FIG. 1. World-line configuration describing quantum AFM 
with the hole; solid (dotted) lines correspond to spin up 
(down) states and the bold fine describes the hole. Arrows 
indicate how periodic boundary conditions are used. 
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FIG. 2. Elementary Monte Carlo updates which form an 
ergodic set. We show updates for and c| end points only 
since procedures for s~ and cj, Cs are identical up to a change 
of notations (i.e., using proper incoming and outgoing lines). 

t tj2 t^J^ t2j3 



FIG. 3. The lowest order (in t and J) identical transfor- 
mation leading to the sign problem: Corresponding graphical 
element changes the sign of the configuration. 

line, which can be traced out and interpreted as spin-1 
magnon excitation. We may now construct an operator 
which has finite overlap with the quasiparticle excita- 
tion as a product s~C|, where s~ is added to cancel the 
extra magnon in the bulk. However, up to a single hop- 
ping transition the above composite operator is identical 
to C|, and we conclude that holes are good quasiparti- 
cles in contradiction with original assumption (probably 
rephrasing the proof of Ref. 0). 
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However, in the presence of the sign problem the above 
consideration should be taken with extreme caution. It 
may turn out that ordered world-line configurations com- 
pensate each other completely and single-configuration 
conclusions are misleading, as suggested in Ref. ^ where 
hole-related sign problem is called the "irreparable phase 
string effect" and argued to cause spin-charge separation. 

The sign problem implies that we may not calculate 
G(t) reliably over long time scales and have to restrict 
our simulation to < 3 4 to suppress sign fluctuations 
by larger statistics. Fortunately, on this time scale GicIt) 
is already in its asymptotic regime and the data are suf- 
ficiently accurate to reveal the ground state properties. 
Formally, calculations are done at finite T but its value is 
more than an order of magnitude smaller than the energy 
of the lowest magnon state in a given system size. For 
each value of J the calculation time was about 2 weeks 
on a PIII-600 workstation. 

In Fig. 4 we show simulated Gko(''') for J/t = 0.4 
and the asymptotic law Z\s_f^e~^^a'^ with the quasipar- 
ticle weight and ground state energy obtained from the 
weight and position of the (5-peak in A\^^{lo). Note, that 
for small values of J the data have to be very accurate 
to describe correctly how G{t) approaches its asymptotic 
behavior Z]^ge~^^o'^ . Error bars are shown but 

are smaller than the symbol size (the relative accuracy is 
better than 10~^ even for points with the largest r where 
the sign problem was the most severe). 

The spectral analysis of Gko(''') was done using 
stochastic optimization procedure developed earlier for 
the polaron problem. A{uj) = J2^=i ^ii'^) is ob- 
tained as an average over spectral densities which opti- 
mize deviations between G(r) and J du!e~'^'^ Ai{uj) . The 
parameter space of Ai{u}) is defined by the step- wise 
constant functions, which, in particular, includes in- 
finitely sharp peaks and is. not associated with any pre- 
defined set of frequenciescil (it is known that J-peaks can 
not be .han dled satisfactorily by the maximum entropy 
methodeea). 

In Fig. 5 we show our results for Ai^g (uj) calculated at 
points J/t = 0.4, 0.2, 0.1. We clearly see a (5-sharp peak 
at the lower edge of the spectrum. The structure of this 
peak is incompatible with the power law singularity since 
its width is smaller than the lowest magnon excitation in 
our system [for J/t — 0.4 the quasiparticle peak width is 
only 0.01 t (!) while the natural scale for the power law 
is set by J] . This is the central result of our paper which 
conclusively rules out spin-charge separation scenario for 
the single hole dynamics in the t — J model and confirms 
finite quasiparticle weight in the thermodynamic limit. 
To verify that finite-size and finite-temperature correc- 
tions are negligible we performed long time simulations 
in a 32x32 lattice at temperature T = J/40 for J/t = 0.2, 
but within the error bars G(r) was indistinguishable from 
the result obtained for L = 16 and T = J/20. 

Unfortunately, the ill-defined problem of numeric ana- 
lytic continuation does not allow to study fine structures 



in the spectral density, especially if they are "screened' 
by low- and high-frequency peaks. (The low-frequency 
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FIG. 4. Gkg{T) (circles) and the asymptotic line 
Zkoe-'^^o^ (dashed) for J/t=0.4. 
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FIG. 5. Spectral functions for J/t = 0.4, 0.2 and 0.1. Fre- 
quency is measured in units of t and the integral J dwA{Lj) 
is normalized to unity. These spectra were obtained for the 
16 X 16 lattice at T = J/20. 

peak is fixed by the asymptotic long-time behavior of 
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G{t) , while the high-frequency peak is fixed by the short- 
time decay of G{t) where the data are extremely accu- 
rate.) Our tests show that multiple peaks in the middle 
can not be resolved by spectral analysis even when we use 
analytically exact G(r) data. It means that the absence 
of multiple string resonances above the ground state in 
our results for A(uj) may not be considered as a proof 
that string potential picture fails in quantum case. We 
would rather consider the second peak as a "course grain" 
description of spectral density at intermediate energies. 
However, if string excitations do exist, their combined 
effect should be seen as the J^/^ scaling law for the peak 
position. In Fig. 6 we plot peak positions as functions of 
(J/i)2/3 for 0.1 < J/t < 1.2 with error bars obtained as 
peak half-widths. We conclude that for the second peak 
the scaling law is obeyed within the error bars. The high- 
energy peak stays roughly at a constant distance from the 
ground state, and clearly the physics behind it is differ- 
ent. 
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FIG. 6. Peak positions as functions of (J/t)^/^. Data 
points for J > 0.4 were taken from Ref. [15] (for J — 0.4 t 
the second peak was not resolved in Ref. [15] because of large 
error bars in G(r)). The two lines are fits y{x) = a + b{J/t)'^^^ 
with a = —3.5 t, b — 3.77 t for the ground state, and b — 5.5 t 
for the first peak in continuum. 
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